NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



MECHANICAL ENGINEERING DEPARTMENT 


(NASA-Ci.<-163026) 'iURBINE BNEViALA 
SW0-C2LIHDiiE PECGEAH ScBlxaiiimdl Report, 1 
Jan. - 1 Jui. VJBO (Connecticut Univ,) 49 g 
HC A03/ai A01 CSCL 21 A 

GJ/37 


li80-2b 6 b5 


Uncids 

27885 



SCHOOL OF ENGINEERING 
THE UNIVERSITY OF CONNECTICUT 

STORRS, CONNECTICUT 


• vr irsairri 


TURBINE ENDWALL TWO-CYLINDER PROGRAM 


Semi-Annual Status Report 
Grant No. NSG 3238 
January 1, 1980 - July 1, 1980 


Submitted to: 

Lewis Research Center 

National Aeronautics and Space Administration 
21000 Brookpark Road 
Cleveland, Ohio 44135 


Principal Investigator: 

Lee S. Langston 
Associate Professor 
Department of Mechanical Engineering 
University of Connecticut 
Storrs, Connecticut 06268 


Introduction 


The following is a report or progress made during the third six-month 
period (January 1, 1980 to July 1, 1980) under NASA Grant No. NSG 3238, 
"Turbine Endwall Two-Cylinder Program." Under this grant an analysis and 
a series of experiments are being carried out to study the three-dimensional 
separation of fluid flow around two isolated cylinders mounted on an endwall. 

The work reported comes under Task I and Task II of the program, and 
deals with: 

a) water tunnel testing 

b) wind tunnel construction and assembly 

c) wind tunnel testing 

d) model design and construction 

e) analytical work 

f) overview of program. 

Water Tunnel Testing 

Water tunnel testing began upon receipt and installation of the turbu- 
lence manipulators discussed in the last progress report^ 

I. Bulk Flow Properties 

Using the hydrogen bubble generator (which is discussed in detail in 
Appendix A) the test section flow was checked for uniformity and a visual 
measurement of turbulence level was obtained. It was found that immediately 
after the perforated plate and screening was installed, the test section 
flow was steady and uniform. However, as testing time accumulated the plate 
and screening became dirty and instabilities began to appear in the test 
section flow. Dealing with this dirt problem requires constant disassembling 
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of the tunnel and cleaning of the turbulence manipulators. Testing at 
low test-section velocities (V < .15 ft/sec) minimizes the instability 
problem. 

II. Main Stream Velocity and Boundary Layer Thickness Measurements 

An unsuccessful attempt was made to measure these low flow velocities 
with a pitot static probe. First, and obviously, the very low dynamic 
heads available for measurement are a problem, and second, we may have run 
into some minimum Reynolds number restrictions associated with our parti cu- 
lar probe, as discussed by Bradshaw . At present, the bulk fluid velocity 
is measured using the hydrogen bubble generator and a strobe operating at 
an equal and known pulsing frequency to obtain a time scale and two fixed 
points in the flow field to obtain a length scale. This method works well 
with velocities as low as V « 0.01 ft/sec. 

One of the best methods available for obtaining a velocity profile 
in the boundary layer is the use of a hot-wire anemometer. We have cali- 
brated a hot-wire anemometer system in the tunnel free stream using the 
previously discussed method for velocity measurement. Using this calibra- 
tion we take our hot-wire probe into the boundary layer to obtain a velocity 
profile. However, the hot-wire probe used is not temperature compensated 
and therefore drifts slowly as the water changes temperature. At present, 
for each velocity profile measured, an independent free stream velocity 
measurement is made using the hydrogen bubble method. We then offset the 
velocity profile given by the anemometer by an amount equal to the difference 
between its free stream velocity measurement and the bubble measurement. 

For this to be valid the anemometer output variation must be linear with 
temperature. 
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It was very quickly realized that a good definition of boundary 
layer thickness is hard to make when interpreting data with some scatter 
(see Figure 1) as well as data where the profile shape changes with test 
section velocity. Clearly, boundary layer thickness alone is not an 
adequate description of what is happening in our boundary layer. More 
information is given by a shape factor with main stream velocity calibra- 
tion. It was found that the boundary layer changes as the upstream turbulence 
manipulators become dirty. In order to categorize boundary layer properties 
we have found it necessary to supply a velocity profile along with each 
set of data taken. Figure 1 shows velocity profiles at three different 
main stream velocities. Over the range of velocities 0.01 - 0.2 ft/sec., 
we have the following ranges of boundary layer thickness, displacement 
thickness, momentum thickness, and shape factor: 

6:. 9-1. 6 in. 

6* : .011 - 0.051 in. 

0 : .007 - 0.012 in. 

H : 1.64 - 3.93 

III. Saddle Point Data 

Figure 2 shows a plot of saddle point location using the coordinate 
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system given in the fir?t progress report . This data was all taken at 
a Reynolds number based on cylinder diameter of Re^ = 1203 while the 
cylinder spacing was first decreased (closing cylinders) and then increased 
(opening cylinders). Dye was injected onto the test section floor through 
both a probe and the dye injection sites in the floor. A Polaroid picture 
was taken of the limiting streamline pattern and the saddle point locations 
were taken from the picture using the dye injection sites as a coordinate 





Figure 1. Velocity Profiles At Cylinder Leading Edge 
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Figure 2a. Radial Displacement of Saddle Point Location as a Function of 
Cylinder Spacing 

Rep = 1203 
D - 1.75 in. 

V = .1163 ft/sec 

6 = 2.0 in. 


Figure 2b. Angula.r Displacement of Saddle Point Location as a Function of 
Cylinder Spacing 
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system on the floor. For some cases the limiting streamlines were "painted 
out" twice and pictures were taken of both to obtain a measure of accuracy 
and repeatability. Then after closing the cylinders they were opened up 
and data was taken that yields a measure of the repeatability of the 
position of the saddle points for a given cylinder separation. 

Figure 2a shows the radial displacement of the saddle point location 
as a function of cylinder spacing. This is similar to the data presented 
in the first progress report^ that was obtained using an open channel 
water table. 

Figure 2b shows the angular displacement of the saddle point as a 
function of cylinder spacing. Again, the trend exhibited by this data is 
similar to the open channel data given in the first progress report. 

There is considerable scatter in the data shown in Figure 2, and as 
such, it should be treated as preliminary. This scatter is due to a 
number of factors: 

1) At lower Reynolds numbers (such as is the case in Figure 2). 
the vortex structure was found to be complex and oscillatory. 

As many as four or five vortices could be observed in the flow 
immediately upstream of each cylinder. Similar flow phenomena 
have recently been reported on by Baker^ for low-speed air flow 
around a single cylinder. 

2) Later testing showed that air could be trapped in the honeycomb 
and screens, leading to a flow distortion in the test section. 

This may have contributed to the scatter shown in Figure 2. 

Consequently, in order to run at higher Reynolds numbers and eliminate 
trapped air, a new inlet section has been constructed, so that before each 
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run the honeycomb and screens can be put in place under water so that 
there is no trapped air bubbles. Tests using this new inlet section will 
be conducted during the next reporting period. 

Mind Tunnel Construction and Assembly 

During the reporting period the construction and assembly of the wind 
tunnel was completed. An overall view of the tunnel is shown in Figure 3. 

Reference kiel and pitot probes were mounted upstream of the test 
section and were connected to a multi tube water manometer, along v/ith 
the static taps on the inlet. The blower motor was wived into the labora- 
tory electrical system and the tunnel was first run on February 12, 1980. 

As was noted in the last progress report, the Flexcore honeycomb for 
the inlet was a long lead-time item, and was not received until April 1980, 
The initial checkout of the tunnel (without the Flexcore) showed that at 
a blower vortex valve setting of about 70®-80° (90° is wide open) the test 
section velocity was 100 ft/sec, just about on design. Thus, testing will 
be possible at Reynolds numbers that are realistic, as was planned. 

After the Flexcore was received it was "sewn" onto the outside of the 
porous plate of the inlet with thin transformer wire. Holes were drilled 
in the top endwall of the tunnel at the middle of the test section and at 
the exit of the inlet, every 6" across the 6 foot span of the tunnel. A 
probe actuator could be mounted at each one of the holes so that a pressure 
and velocity profile could be made, endwall to endwall. In addition, a 
set of stairs and a platform was constructed on top of the test section 
to provide ready access to the probe actuator. 

Since the tunnel has been in operation, it has been found that the 
blower vibrates the test section to some degree. A rubber diaphragm had 
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Figure 3. Overall view of wind tunnel showing honeycomb mounted on porous plate at inlet, Stanitz inlet nozzle 
strainht duct, test section with lad**er. dump diffuser and blower. 
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been used to connect the blower to the diffuser, but this has not 
proven to be up to the task of eliminating all vibration. Some vibration 
is also transmitted through the floor. 

In the last few weeks, steps have been taken to eliminate the blower 
vibration by mounting the blower on springs and installing a flexible 
duct between the blower and the diffuser. To a first approximation the 
transmissibility , T, of the spring-blower-mass system is given by 


T 





0 ) 


where damping has been neglected and where 

T = transmissibility, i.e. the force due to the blower vibration 
transmitted through the springs to the floor, divided by 
the same force for a rigid mounting. 

( 0 ^ - the frequency of the blower vibration (roughly 1800 rpm). 

Wg - the natural undamped frequency of the blower mass-spring mount 
system. 

The last term, «g is given by 



where 

K - spring constant 
- blower mass 


It can be seen from (1) and (2) that to minimize T, the mass of the blower, 
mg, should be large. One can do the same thing with a smaller value of 
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K, but there is a practical lower limit on K. [If the springs are too 
soft (low K) the blower will have a tendency to "jump" around]. Thus, to 
lov:(sr T, about 3000 pounds of concrete were poured into the base of the 
blower. When this cures (about mid-July 1980), the spring mounts and 
the flexible duct will be installed, and the concrete base and blower 
will be raised off the floor. 

Wind Tunnel Testing 

During the reporting period an extensive test program was undertaken 
to check out the characteristics of the wind tunnel, before any model 
testing was done. 

First the static pressure distribution of the Stanitz inlet was 
measured, using the wall pressure taps that had been installed for that 
purpose. All pressures were measured relative to the static pressure of 
the pitot static reference probe upstream of the test section. The 
pressure differences were read using a Flow Corporation micrornanometer 
that is accurate to + 0.0002 in. of manometer fluid (N-butyl alcohol). 

Figure 4 shows the pressure distribution that was measured along the 
centerline of the inlet, as a plot of q/q^ (local velocity head, 
over reference velocity head, ) vs. surface distance along the inlet. 
The reference velocity was calculated from the reference total pressure 
measured in the inviscid portion of the flow upstream of the test section 
and the ‘static pressure measured at the point x = 0.5 ft. at the lower 
surface of the inlet. This value was used rather than the velocity measured 
at the test section because the latter is slightly higher due to the dis- 
placement thickness effects of the tunnel boundary layers. 
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Figure 4. Acceleration Through Nozzle - Comparing Stanitz 
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Figure 5. Nozzle Cross-Section — 6 in. Upstream of Nozzle Exit Plane 



Also plotted in Figuna 4 is the Stanitz potential flow pressure 
distribution that the inlet was designed from. As can be seen, there 
is good agreement with potential flow on the nozzle inlet centerline for both 
the top and bottom of the inlet. 

Figure 5 shows the static pressure distribution across the exit of 
the nozzle inlet on the bottom surface. The measured values of static 
pressure show some variation (note that ordinate scale has been expanded) 
away from the tunnel centerline. While this variation is not large, its 
presence is disturbing, since the inlet was carefully designed and con- 
structed to provide uniform flow. More will be said about this variation 
in what follows. 

Total pressure profiles were measured midway through the test section 
and upstream, at the exit of the inlet, using a United Sensor -g" head kiel 
probe, with an acceptance angle range of +45°. All pressure differences 
were measured with the micromanometer described above. The probe was 
positioned using an L. C. Smith probe actuactor. Model No. BBR-123-360. 

The actuator had been calibrated for position and yaw angle before the 
tests. Using the digit digital panel meter on the probe actuator 
indicator, the probe position could be read to + 0.02" and the probe yaw 
angle to within + 0.4°. 

The test section has a 1x6 ft. cross-section. At the test sec- 
tion,. traverses of total pressure were taken from floor- to-ceiling (the 
1 ft. dimension) at locations 6 in. apart from the north side to the south 
side of the tunnel (the 6 ft. dimension). 

Three of these total pressure traverses are shown in Figure 6, where 
distance from the tunnel floor is plotted against the total pressure 
coefficient, Cp^, given by 
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Figure 6, Total pressure profiles at the test section. 
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Cp. = 


Pt - 




where - measured total pressure 

Pf^ - reference total pressure (always measured with reference 
probe on north side of tunnel). 


The data in Figure 6 shows that the in viscid core of the flai; in the 
test section is fairly uniform in total pressure except at the south side 
of the tunnel, where the total pressure is slightly higher (about 1-3% 
which amounts to approximately 0.01-0.03 in. of water) than the centerline 
and north side total pressures. Secondly, the tunnel floor boundary layer 
gets thicker from the centerline to either the north or south, while the 
ceiling boundary layer is thinnest on the north wall. Thirdly, the data 
in Figure 3 shows that the ceiling boundary layer is larger than the floor 
boundary layer. 

The slightly non-uniform static pressure distribution at the inlet 
exit, the slightly non-uniform total pressure distribution at the test 
section and the asymmetric boundary layers on the test section floor end 
ceiling are not to be expected, considering the care that was taken In the 
design and construction of the inlet. 

The internal surfaces of the tunnel were inspected and it was found 
that the tunnel ceiling had one joint that caused a small front facing step. 
This was corrected ‘-ut no change in the flow field was noted. 

In an attempt to make the floor and ceiling boundary layers more alike, 
a trip wire was installed at the exit of the inlet. The wire was D = 0.02 
inches in diameter and was chosen based on a trip Reynolds number of 

UcoD 

~ 900 
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based on the criteria given in Schlichting . Total pressure measures were 
again taken at the test section and showed that the trip wire had little 
or no effect on either the ceiling or the floor boundary layers. Another 
test was done using a larger trip wire (D = 0.125") which also showed no 
effect on the test section boundary layers. 

Traverse holes ware !:han cut at the exit of the inlet, and traverses 
were taken every 6 inches (in the 6 foot dimension) with both the keil probe 
and a United Sensor pitot probe. Data from three of these traverses is shown 
in Figure 7. Figure 7a shows velocity and Figure 7b shows total pressure. 

The velocity plot shows that there is a slightly higher velocity in the 
south side of the tunnel. The total pressure traverses show essentially 
the same features as test section traverses, but scaled down, since the 
inlet traverses are about 6 feet upstream of the test section traverses. 

Thus far, the conclusion is that the flow distortion and boundary layer 
asymmetry is coming from the inlet. Efforts are now under way to see if 
these effects are (a) being generated by the inlet itself (porous plate 
honeycomb, etc.), or (b) entering the inlet from the room itself and are 
not being sufficiently smoothed out by the porous plate. 

Model Design and Construction 

During the reporting period the first pair of cylinder models for the 
wind tunnel were designed. Construction of these models is now in the 
final stages. A photograph of one of the cylinder models is shown in Figure 
8. The model is an aluminum cylinder that was accurately machined on a 
lathe to an outer diameter of 6.25 inches. 

The two cylinders will be placed side-by-side in the wind tunnel test 
section by inserting each of them through an airtight fixture in the top 
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Wind tunnel cylinder model. At left is 6.25" OD aluminum 
cylinder with 59 static pressure taps, and mounting fixture 
at <^ioht is pressure tao manifold. 
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of the test section. A clamping device (see Figure 8) will hold them in 
place and will press each cylinder against the floor of the test section. 

An 0-ring at each cylinder base will seal the cylinder - test section 
floor junction. 

Each cylinder has 59 pressure taps, with the greatest number of taps 
near the base of each cylinder, where the separation process will occur. 

The cylinders and the clamping fixture are constructed in such a way that 
a cylinder can be rotated (indexed) to change the relative circumferential 
location of the 59 pressure taps. 

The pressure taps are connected on the inside of the cylinder to a 
manifold of 0.047" ID stainless steel capillary tubes. The use of this 
manifold (shown in Figure 8) allows the model to be constructed from a 
single, more sturdy, cylinder, rather than cutting the cylinder in half 
and individually soldering each capillary tube to each pressure tap. 

The spacing between the cylinders can be varied by use of various 
size plexiglass spacers in the top of the test section. Since the bottom 
of each cylinder has a compressed 0-ring seal, the floor of the test 
section can be kept clean (no junctions or seams) for the flow visualiza- 
tion experiments. 

The construction of these first cylinder models should be completed 
by mid-July 1980, at which time flow visualization and pressure measure- 
ments will be done for a range of cylinder spacings. 

Analytical Work 

A two-dimensional potential flow program obtained from Professor Patrick 
Kavanaugh of Iowa State University was adapted to the UConn Computer System. 
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This program is based on the distributed source-sink method of Smith and 
Pierce. It is an on-body program, i.e. as output, it provides pressure 
distributions on the body itself. The output of this program will be used to 
compare with the experimental data from the pressure taps on the cylinders, 
to show any departure from potential flow. 

Off-body pressure will also be taken and it is necessary to compare 
this data with a potential flow calculation. Thus during the reporting 
period, an off-body two-dimensional potential flow program for flow around 
two cylinders was completed. A discussion of this program is given in 
Appendix B. This program will be used to: 

a) provide an early first-order correlation of saddlepoint 
location (gotten from flow visualization data) and potential 
flow pressure gradients at the saddlepoint location. 

b) compare to flew field experimental data. 

A third computer program has also been worked on and adapted to the 
UConn Computer System. This program was also obtained from Professor 
Kavanaugh, and is a general, two-dimensional contour plotting routine to 
be used in conjunction with a Cal comp Plotter. This program gives us the 
ability to plot isobars, streamlines, etc., be it from experimental data 
or from the output of a potential flow program. 

Overview of Program 

As noted in the last report^ the water tunnel tests are behind schedule 
due to the inlet distortion problems that were encountered. However, the 
tests are far enough along to give us a good feel for cylinder spacing, wall 
effects, and low Reynolds number behavior, all of which have impacted wind 
tunnel model design. 
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The wind tunnel testing is slightly behind in schedule. The cylinder 
models will shortly be ready to be installed in the wind tunnel to 
check them out, but no data will be taken until the wind tunnel flow has 
been made more uniform. 

Once better uniformity is achieved, testing will begin. The first 
tests will be to vary cylinder spacing and get the location of the saddle 
points by using flow visualization techniques. 

The off-body potential flow program will be used to calculate normal 
and streamwise pressure gradients at the saddlepoint locations. This will 
give a first order correlation. Detailed flow measurements will then follow. 
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APPENDIX A 

The Hydrogen Bubble Generator 

Introducti o n 

The purpose of this project was to produce a hydrogen bubble generator 
for use as a flow visualization medium in a flow of water. Herein we 
present an explanation of the chemical processes involved in bubble 
generation as well as the design of a simple device for that purpose. 

Background 

In many fluid iiiechanics investigations there is a great need for 
visualization of comple:, fluid patterns. Flow visualization provides a 
means for exploring some flow phenomena that are not predictable by mathe- 
matical theory. Using hydrogen bubbles as tracer particles in a water 
flow is a good visualization technique. The first to apply this technique 
Was a German scientist named F. X. Wortmann^ who was followed in the 
United States by E. W. Geller^*^. Clutter and Smith^ of Douglas Aircraft 
Company completed extensive research perfecting the hydrogen bubble 
visualization technique. 

Theory 

The hydrogu.l? bubbles are produced by the electrolysis of water and 
must be of a size and quantity so as not to disturb the flow field. The 
electrolysis of water is a simple reaction which decomposes water into its 
basic elements; 

HgO *■ Hg + ^0^ 

This reaction does not occur spontaneously therefore energy must te supplied. 
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As seen in the chemical equation above, there are twice as many moles 
of hydrogen produced than there are moles of o>^gen. By applying 
Avogadro's Law at constant temperature and pressure it can be seen that, 
volumetrically, there are twice as many hydrogen bubbles produced as 
there are oxygen bubbles. It is for this reason that we choose to use 
hydrogen bubbles rather than oxygen bubbles to trace the flow. 

The reaction occurs as follows; for the simple electrolysis cell 
shown in Figure A-1, bxygen bubbles are produced at the positive terminal, 
the anode, while hydrogen bubbles are produced at the negative terminal, 
the cathode. 

The choice of a material (metal) to be used for terminals is based 
upon the chemical reactions that take place at the terminals as well as 
a strength requirement. (The Hg source or cathode must withstand drag 
force of fluid). For instance, from the chemical reaction standpoint, 
platinum is a very good material. 

The two reactions which take place at a platinum cathode (See Figure 
A-2) are; 

—1 ^"2 AGi 

+ 4e + Pt^ 5 j >> - P*^2(S) 

and 


2H 


,+i 

(L) 


+ 2e 


-i 




Note that the first reaction yields an undesirable solid precipitate, 
platinum hydrate, while the second yields the desired hydrogen gas. Both 
reactions require energy for them to proceed to the right. The energy 
required is called the free energy of formation (AG) and is much less for 
the gaseous hydrogen reaction than it is for the platinum hydrate reaction, i.e. 


AG2 <« AG^ 
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Both reactions will occur but for all practical purposes no platinum 
hydrate will bo formed. Platinum is also very corrosion resistive. 

However, in spite of all these advantages, platinum is limited In its 
use by its weakness. Other metals such as tungsten have been found to 
be an acceptable substitute. Tungsten has good corrosion resistance 
properties, along with good strength and a reasonable price. 

At the anode we also have two chemical reactions. Again, platinum 
would make an ideal material for the terminal; however, a less expensive 
metal such as copper can be used (Figure A-3). 

With a copper anode the two possible reactions are: 

I AGo 

°(L) ^^(S) ^ ® 

and 

-2 , “1 

The free energy of formation of these reactions are close in value favoring 
the reaction which forms solid copper oxide. This will produce a large 
decay of the copper anode along with oxygen formation. This is not a 
great disadvantage in that we can probably install the anode in some 
unimportant region of the flow. 

Application 

Now that we have presented a methodology for hydrogen bubble formation 
we can discuss its application to flow in a water table. 

The bubbles are introduced into the flow field by a carefully positioned, 
thin wire which will be designated as the cathode and energized accordingly. 
The anode terminal must then be placed anywhere in the water table in order 
to complete the cell. The anode is usually placed downstream of the test 


26 


section so as not to disturb the test section flow field (See Figure A-4). 

3 

In their research, Clutter and Smith found that at flow velocities 

below 5 feet per second, wire diameters up to .005 inches could be used 

with little flow interference. Figures A-5, A-6, A-7, and A-8 are a 

3 5 

graphic presentation of data taken by Clutter and Smith and Mattingly 
during their flow visualization studies. 

3 

To produce a continuous supply of bubbles. Clutter and Smith found 
that a supply voltage of 25 to 400 volts is necessary, depending upon the 
resistance to current flow associated with the water table (electrolytic 
cell) being used. For our particular table, it was found that a voltage 
between 100 and 150 volts was sufficient to produce a continuous supply 
of bubbles adequate for good flow visualization. This voltage range 
yields a current range from .25 to .5 amps on our water table. The voltage 
and therefore current levels must be adjusted to give a "good" stream of 
bubbles (there must be enough bubbles so that they can be seen well, and 
they must not be so large that they float up) for each flow velocity tested. 
Currents can be increased for a given voltage level by seeding the water 
(table salt). 

By pulsing the supply of voltage to the cell, a cyclic generation of 
hydrogen along the cathode wire produces patterns like those of Figure A-9. 
By measuring the pulse duration electronically a velocity profile can be 
obtained simply by measuring the pulse width. There are several pulse 
generators that can supply this signal, but not at the voltage and current 
required. To obtain the high voltage square wave (Figure A-10) a signal 
generator is used to activate an electronic switch to an additional power 
supply (Figure A-11). (A complete schematic is shown in Figure A-18.) 
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Figure A 



. Flow Through a Converging Channel, taken from r\eference 7. 
Upper Photo shows Time Lines Made by Hydrogen Bubbles from 
Pulsed Wire. Lower Photo shows Combined Time-Streak Markers 
Produced by a Pulsed, Partially Insulated Wire. 
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It is necessary to have a sharp square wave of voltage in order to 
produce a consistent bubble size throughout the pulse width. Without the 
sharp square wave of voltage, the bubble size will vary throughout the 
pulse width causing large separation of the bubbles as the bubble sheet 
propagates through the test section (see Figure A-12). It can be seen 
in the figure that the slight variation in voltage causes inconsistent 
bubble sizes. The buoyancy and drag forces will be different for each 
bubble size causing the larger bubbles to move slower and rise more, while 
the small bubbles follow the flow stream. 

Stoke 's analysis of the forces on a buoyant sphere in a steady low- 


speed water flow, shows a buoyancy to drag ratio of: 

B/D = 




18Uv 


where 


d = bubble diameter 
U -• flow velocity 
V = kinematic viscosity 
B = buoyancy force 
D = drag force 


It can be seen in Stoke 's analysis that as the flow velocity decreases, or 
bubble diameter is increased, the buoyancy to drag ratio will increase, 
causing the bubbles to deviate from the flow field. A free body diagram on 
an individual bubble is given in Figure A-13. Bubble diameters between 
0.0005 and 0.001 inches will have a rise rate due to buoyancy of about .05 
inches per second. With a flow velocity between 0.1 and 5 feet per second 
this rise rate is insignificant. 

Flow velocity and the diameter of the cathode wire play another role 
in determining bubble size. The faster the flow, the greater the drag forces 
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Figure A-12. At top, an edge view of pulsed time lines showing the effect of 
uneven bubble sizes. At middle and bottom, a top view of pulsed 
time lines showing a non-square wave behavior. Taken from Reference 8. 







Figure A-13. To accurately follow the flow the buoyance to drag ratio must be small. 
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on the hydrogen gas around the cathode wire. This drag causes the gas 
to separate quickly from the wire thereby producing small bubbles. 

Schraub and Kline®, at Stanford University, perfected the techniques 
necessary to detect path lines, streaklines, and streamlines in low speed 
water flows (Figures A-14, A-15, and A-16). 

Figure A-17 shows some possible configurations for lighting the stream 
of bubbles. It has been found that the 125® angle between the light source 
and the line of vision and the flat black color used for contrast, are 
important in getting good photos of the bubbles. Figure A-18 shows the 
circuit diagram of the electronics for the bubble generator, and Figure 
A-19 shows the bubble generator control panel. 
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Figure A-19. Bubble Generator Control Unit and Power Supply. 
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APPENDIX B 

Two-Dimensional Off-Body Flow Program 


Summary 

The two-dimensional off-body potential flow program was developed 
to provide velocities* pressures, and streamlines in the flow field 
around arbitrary two-dimensional bodies. The bodies used in this program 
may be any series or arrangement of closed bodies, with lift and/or in 
an infinite cascade. The results are printed in a table for various 
points in the flow field. Options also exist for plotting isobars or 
streamlines on the Cal comp Plotter. Fundamental to the two-dimensional 
off-body potential flow program is the on-body potential flow program, 
as the parameters needed for the off-body program are generated by the 
on-body program. In addition to these two programs, a third program 
used in making plots of isobars and streamlines is a contour plotting 
program. This program and the on-body program are products of previous 
work at the University of Iowa. 
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Symbols 


u 

X velocity 

V 

y velocity 

'p 

stream function 


source or sink of strength 

U 

unit free stream flow 


strength of circulation 

npts 

total number of all body points 

nobods 

number of bodies in the flow field 

s 

length of line singularity 

(Xp.yp) 

arbitrary point in the flow field 

("p*yp) 

translated point in field 


midpoint of line source segment 

0 

angle of line source to horizontal 

s 

pressure coefficient 

sp 

cascade spacing 

Re 

r«!al part of a complex function 

Im 

imaginary part of a complex function 

Zi,Z2 

complex vectors 
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Theory 

The basis for potential flow calculations is the solution to Laplace’s 
equation. Laplace's equation is derived through the irrotationality 
condi ti on ; 

au. ly. - n 
ay "ax " ^ 

In addition the stream function is defined as: 


u = M 
3y 


= 

" ax 


which gives after substitution into the irrotationality condition; 


L- (ilk) _ Jl. (_ ii) = 0 
ay '■ay' ax '■ ax' 


^ = 0 


ay' 


ax' 


This is Laplace's equation. A solution or set of solutions to Laplace’s 
equation can be made into a sum which is also a solution. This condition 
of linearity facilitates the calculation technique used in this program. 

Bodies in the flow field are represented by distinct points describing 
their contours. The program uses these points to generate line source- 
sink singularities which adjusts the free stream flow to that which 
satisfies the Neumann condition on the boundary of the bodies. In the 
case of lift the flow will satisfy the Kutta condition, also. The overall 
approach here is based on the work of Hess and Smith\ 

The stream function for a single source at the origin is given by: 

ip - 0^ tan"V/x 


42 


Adding a free stream flow yields 

i|) = Uy 

and the contribution from circulation 


^ ^ lo9g(x^+ y^)^ 

gives the stream function for a single source located at the origin, 
with lift in uniform flow. 


ip = Uy + a^tan"V/x - ^ ln(x^+ y^)^ 

For every body point specified there is a corresponding strength 

of the source-sink singularity, To account for all the body points 

a summation may be taken. 

npts p , , 1.. 

= Uy + I (cr.tan y/x - ^ ln(x + y )'^) 
i=l ’ 


Although the circulation value is similar for all points in the same body, 
if there is more than one body the circulation must be summed, yielding: 


npts 


nobods 


= Uy + I (o.tan"V/x - I ^ ln(x^+ y^)^) 
i=l ^ K=1 


The above value of the stream function is true for all different 
singularities, all located at the origin, however. Therefore each source- 
sink must take on the value of a line singularity, removed from the origin, 
and angled such that the sum total of all singularities maps out the 
bodies in the flow field. To develop the line singularities consider 
an element at the origin orientated as shown. 
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The stream function for the element can be found by integrating the 
source-sink function over the region -s/2 to s/2. 

,s/2 


ip = [ a.tan~ 
J-s/2 ^ 


dx 


Making the substitution: x' = Xp-x, dx' = -dx, when x - -s/2, x' = Xp+s/2 

and when x = s/2, x' = Xp-s/2. The stream function becomes: 


/p-s/2 

ip - a. tan -?■ dx' 

Jx+s/2 ^ ^ 


Xp-s/2 _ , 

= - [-a^tan ^ ~+ o.7r/2]dx' 

Jxp+s/2 ’ ^ 


Since the value for is constant over each element, it may be taken 
outside the integral. 




Xp+s/2 

= a. f [tt/2 - tan“* ^ ] dx' 

^ JXp -s/2 
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1 Vs/2 , Xn"S/2 

tp = a^U/Z s -(Xp+s/2)tan y " " •♦• (Xp-s/2)tan — + 


r + (Xp-S/2)= ^ 


A coordinate transformation must also be used, as each segment is 
not located at the origin with zero angle of orientation. The segment 
shown gives an example of how this transformation is accomplished. From 
the diagram it can be seen that subtracting x^^ and from the midpoint 
will bring the line segment center to the origin. A rotation, to the 
horizontal of the line segment may be accomplished by an angle, -0. The 
angle, e, is found from geometry and is equal to 

0 = tan"^ ~ ^ni 
Xj - Xpi 



The correcting displacements to Xp and yp are achieved by: 
Ay = s/2 sin 0 Ax = s/2 cos 0 
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The translation and rotation are applied to the point (Xp, and 


(x^,yp); with: 


- 


x_ - X„ “ Ax 


ppm 


y;, = yp - ^m - 

Adding the effects due to free stream and circulation, the stream function 
now becomes 

npts 1 x*+s/2 V -1 

= ^yn^l Co.(tt/2 s -(x'+s/2)tan -P - r - ^ (x* -s/2)tan -Pp— + 

P i=l ^ P P 


Yr. +(x'+s/2)^ nobods u 

^ +(x;,-s/2)" K-1 P P 

Now to calculate the velocities of the off-body points, apply the 
formulas 

P ■ ly P"'' P = - • 


npts 1 (x!+s/2)^+ y'^ nobods 

^ I (cr. ^ In ”T .2 - I 

i=l ' 2 (x!-s/2)^ K=1 


u = y. 


2it 


, .1 x;,+s/2 _J -s/2 

(tan tan )) 


npts ,x;,+s/2 x;-s/2 nobods (x;+s/2)='+y^" 

? rx-„-l P 4-r.^, ^ _E 1 4- V In K V. 


V = 


(a^. [tan' 


•i_E 


-i — - tan 


I 


'ii In 


K=1 ^ (Xp-s/2)^+y^^ 


)) 


i=l ' -"P "P 

These are the formulas used in the calculation for velocity and 
pressure in the program. C , the pressure coefficient is calculated from 


r - 1 v')' 
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The velocities for a point in a cascade flow are calculated somewhat 
differently than the above formulas for uniform lifting flow. Calculations 
are arrived at in the complex plane. In the report for the on-body program, 
the author derives the velocity formulas using a complex potential function. 
For uniform flow and lift the formulas are given below, where sp equals 
the specified cascade spacing and a is the angle of the line source element 
with the horizontal. 


u = 




npts 

L V / 


npts 
V = I (-o.Im[e 
i=l ^ 


-la 


ln( 


sinh 

sinh 


TfZi 


J£_) 

TTZz ' 

sp 


nobods 

]+ I 

K=1 


Ztt 


Re[e"'“ln( 


sinh 

sinh 


TTZg 

sp 


■) ]) 


Zj = (Xp - *i) + i(J'p - yi) 

^ ' (’‘p - >< 2 > * '■% • 

It remains only to calculate 0^ and to complete the analysis for 
potential flow. The purpose of using the on-body program with this program 
is to generate and as the on-body program uses these values, also. 

The derivations and formulas for arriving at 0^ and Fj^ are spelled out 
in the on-body report and are far too numerous and complicated for the 
context of this report. Due to the use of complex numbers for cascade flow, 
it was not as easy to calculate the stream function analytically as with 
uniform flow. Therefore, unfortunately, there is no provision for plotting 
streamlines for cascade flow. The streamlines one obtains are those simply 
for the repeated configuration entered into the program, taken as if it 
were alone in the flow field. 
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